global alpha beta;

alpha = 0.0544;
beta = 0.0106;

x = 0:100:20000;

K = [-2500000,0,2500000];

plotConfs = ['k: '; 'k- '; 'k--'];

hold on;

for i=1:3;
    y = real(sqrt(((alpha * x.*x)-K(i))/beta));
    plot(x, y, plotConfs(i, :));
end

xlabel('x');
ylabel('y');
legend('K < 0', 'K = 0', 'K > 0', 'Location', 'NorthWest');

hold off;
